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Revisiting the angular momentum growth of protostructures 
evolved from non-Gaussian initial conditions 
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ABSTRACT 

I adopt a formalism previously developed by Catelan and Theuns (CT) in order to estimate 
the impact of primordial non-Gaussianity on the quasi-linear spin growth of cold dark matter 
protostructures. A variety of bispectrum shapes are considered, spanning the currently most 
popular early Universe models for the occurrence of non-Gaussian density fluctuations. In 
their original work, CT considered several other shapes, and suggested that only for one of 
those does the impact of non-Gaussianity seem to be perturbatively tractable. For that model, 
and on galactic scales, the next-to-linear non-Gaussian contribution to the angular momen- 
tum variance has an upper limit of ~ 10% with respect to the linear one. I find that all the 
new models considered in this work can also be seemingly described via perturbation theory. 
Considering current bounds on /nl for inflationary non-Gaussianity leads to the quasi-linear 
contribution being ~ 10 - 20% of the linear one. This result motivates the systematic study 
of higher-order non-Gaussian corrections, in order to attain a comprehensive picture of how 
structure gravitational dynamics descends from the physics of the primordial Universe. 
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1 INTRODUCTION 

Gravitational instability and hierarchical growth of Cold Dark Mat- 
ter (CDM) density perturbations provide an elegant description for 
the acquisition of angular momentum by protostructures. Accord- 
ingly, patches of matter are spun up by tidal t orques exerted by the 
surrounding Large-Scale Structure (LSS, se e Peebles 1969 Jl97l|; 



lDoroshkevichlll97q lWhitdll984l: Heavens & PeacocHll988l) . At 
the linear level the spin growth is given by the coupling of the first- 
order deformation tensor and the inertia tensor of the patch, and 
it agrees q uite well wit h more accurate numerical results (see the 
review bv (Schafej2009l and the references therein). 

Acquisition of angular momentum beyond the simple linear 
description has been tackled by means of Lagrangian perturba- 
tion theory in a series of paper s by Ca telan and Theuns dCatelarj 
ll995l:ICatelan & TheunJl996j |rl [T997l . CT henceforth) almost two 
decades ago. In particular, they found that the next-to-linear cor- 
rection to the growth of ensemble-averaged spin is non- vanishing 
only if primordial density fluctuations are non-Gaussian. CT ex- 
plored several non-Gaussian models, and concluded that only for 
one of these does the angular momentum acquisition appear per- 
turbative ly tractable: the log -normal model for the gravitational po- 
tential of Moscardini et al. For this template, and adopting 
a representative mass scale of M ~ 10 12 /r' M G (with h = 0.5 and 
Gaussian filtering), CT found an upper limit of ~ 24% for the quasi- 
linear non-Gaussian contribution to the spin variance. This figure 
translates to a ~ 10% value when rescaled to a smaller and more 
typical galactic mass of M = W w h- l M s (with h = 0.7 and a real- 



space top-hat filter). For other templates the non-Gaussian contri- 
bution is comparable to, or larger than, the linear term, suggesting 
the impossibility of a perturbative expansion. 

Since then, the issue of angular momentum growth in non- 
Gaussian cosmologies has not been investigated further. On 
the contrary, new and more general models of primordial non- 
Gaussianity exist nowadays and, most importantly, constraints on 
the level of primordial non-Gaussianity coming from the Cosmic 
Microwave Background (CMB) and the LSS have dramatically im- 
proved over the last decade. Gi ven the cosmological relevance of 
primordial non-Gaussianity (see Barto lo et al.l2004llChenl2010l for 
recent reviews) and the significance of the CDM halo angular mo- 
mentum acquisition for the formation and evolution of galaxies, it 
is important to update this topic. Specifically, it is interesting to 
explore the amplitudes and behaviors of the quasi-linear contribu- 
tions to spin growth given by non-Gaussian models that are popular 
nowadays. This is the scope of the present letter. 



The rest of the manuscript is organized as follows. In Sec- 
tion [2] I review the linear and next-to-linear contributions to the 
ensemble-averaged spin growth of matter patches. In Section [3] I 
summarize the non-Gaussian cosmologies explored here. In Sec- 
tion|4]results are displayed and in Section[5]conclusions are drawn. 
Where needed, I adopted the following cosmological parameters: 
fi m .o = 0.272, fi A , = 1 - n m>0 , fib.o = 0.046, H = 100/i km s~' 
Mpc -1 with h = 0.704, cr 8 = 0.809, and n s = 1. 
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2 ANGULAR MOMENTUM ACQUISITION 
2.1 Lagrangian displacement 

In Lagrangian theory the comoving position x of a mass element 
at time r can be written in terms of its initial position q and a dis- 
placement vector field S, as 



x(q,T) = q+S{q,r) . 



(1) 



Following CT here I used a time var iable r that is re lated to the 
standard cosmic time t by dr = dt/a 2 ( Shandarin 1980), where a is 
the scale factor. 

Perturbative approximations to this exact expression can be 
found by expanding the displacement field in a series, 



(2) 



where Si corresponds to the Zel'dovich approximation Si(q,r) = 
D(j)Vifi\{q). Here D(r) is the growth factor of linear density per- 
turbations, which in a Einstein-de Sitter universe reads D(r) = r~ 2 . 
The fu nction if/j is the first order (Zel'dovich) displacement po- 
tential feel ' dovichll 1 97dh , related to the linear density perturbation 
field by the Poisson equation Aif/ t (q) = S(q), so that in Fourier 
space tjriip) = S(P)/P 2 - 

The second-order term of the displacement field can also be 
separated in time and space, according to £2(9, t) = E{ryj)j/ 2 (q). 
The growth factor E(r) reads E(r) = -3t~ 4 /7 in an Einstein-de 
Sitter universe, while for its more general expression I refer to CT. 
The second-order displacement potential can be related to its first- 
order counterpart in Fourier space by 

hip) = -4 r ^^k^?8i i ip l +p 2 -p)\K{p l ,p 2 )x 

p 1 Jm* (Sir L J 

x faipOfa&i). (3) 

In the previous Equation A"(p, , /?,) is a symmetric integration kernel 
defined as 



K(PuP 2 ) = I \p\p\ - (Pi -Pi?} = \p\p\{} 



(4) 



where yu is the cosine of the angle between the two wavevectors p, 
and p 2 . 



2.2 Spin growth 

The angular momentum of the matter initially contained in a co- 
moving Lagrangian patch T of the Universe at time r can be written 
as an integral over T, 



J(r) = a 3 (T)p m-0 J dq[q+ S(q, r)] X 



dS(q,r) 
dr 



(5) 



By considering the series expansion of the Lagrangian displace- 
ment field S introduced in Eq. l[2}, the angular momentum of the 
patch can be similarly written as 



/ = £/„ 



(6) 



By expanding the Zel'dovich potential around the center of 
mass of the patch (assumed to be, without loss of generality, the 
origin of the reference frame) up to the second order, the previous 
equation takes the compact form 

dD(r) 



dr 



- E aPyO lJ3tT?cry( T ) 



(8) 



In the previous equation £ a p y is the fully antisymmetric Levi-Civita 
tensor, is the Zel'dovich deformation tensor, 



£> 



dq p dq a 



f dp 



PpPa^l(P) ■ 



while la-y is the inertia tensor of the patch, 



Iayif) = « 3 (t)p, 



1,0 J q^q y 



(9) 



(10) 



Summation over repeated indices is implicit. 

Likewise, the second-order term in the series expansion of the 
angular momentum reads 



1 dEM r 
Ji(t) = a s {r)p m , Q —^- j dqqx Vifr 2 (q) , 



(11) 



which, under a second-order Taylor expansion of the displacement 
potential takes the form 



(12) 



It can be shown that the second-order deformation tensor in Fourier 
space reads 



2 JScr 



r dpjd 



2 K (PuPi) x 



(13) 



r) e IIP1+P2 

x fi(Pi)h(Pi) > 

in terms of the Zel'dovich potential. 
2.3 Ensemble averages 

In order to simplify the previous results, it is meaningful to consider 
the ensemble average of the square of the angular momentum. It 
then follows that, up to the next-to-linear order, 



IIJ(r)|| 2 )- ||J 1 (t)|| 2 ) + 2(/ 1 (t)-/2(t)) 



where 



(ll/i(r)|| 2 ) = - 



dD(r) 



dr 



(y\ - 3v 2 ) ( 



(14) 



(15) 



In the previous Equation a M is the mean deviation of the mat- 
ter density field smoothed on a scale corresponding to mass M, 
while V] and v 2 are the first and second invariant of the inertia 
tensor, respectively. To be more precise, if A\, A 2 , and A3 are the 
three eigenvalues of the inertia tensor, then Vj = A\ + A 2 + A3 and 
v 2 = /1|/L + ~t" A 2 Aj. 

The next term is non- vanishing only if density fluctuations are 
non-Gaussian, and reads 



,1^1,^ 2 dD(T) dE(T) ( P- ■> \ 
</i(t) ■ J 2 (t)) = —^z — y\ - 3v 2 )w M • 



15 dr dr 



(16) 



where 



The first-order term of the angular momentum series takes the form 



Ji(t) = a (r)p, 



dD(r) C 



xV^( 9 ). 



(7) 



,< r d Pl d 

u M = -15 -— 

Ju 6 (2tt 



i d Pi 

T )6 



\\p l +P 2 \\ 2 K(p 1 ,p 2 )W 2 R (\\ Pl +p 2 \\)X 



(17) 
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In the previous Equation represents the bispectrum of the 
Zel'dovich potential, which is now explicitly smoothed on a scale 
R = (2GM/Q. m ,oH(,y /3 ■ I assumed the standard real-space top-hat 
smoothing. The Zel'dovich potential can be related to the standard 
gravitational potential <p by making use of the Poisson equation, 

fa(P) = %i£j^ftP)eFWWP). (18) 

where 7"(p) is the cold dark matter transfer function dBardeen et al.l 
1 19861: ISugivamj[T99^) . The integral in Eq. (T7J has to be solved 
numerically for realistic bispectrum shapes. Fortunately, the bis- 
pectrum usually depends only on the magnitude of its three argu- 
ments, so that the above six-dimensional integral reduces to a three- 
dimensional one. 



3 NON-GAUSSIAN SHAPES 

I considered five different shapes for the primordial bispectrum, 
that are briefly described below. The first four are motivated 
by inflationary physics, and the amplitude of non-Gaussianity is 
given by the parameter / NL (assumed to be constant). The fifth 
is non-inflationary in nature, and hence independent on /nl. See 
iFedeli etalJboilh and references therein. 

3.1 Local shape 

This bispectrum shape arises when a light scalar field, addi- 
tional to the inflaton , contri bu tes to the curvature perturbations 
jBernardeau & Uzanl |2002|; iBabich. Creminel li. & Zald arriagal 
|2004 ISasaki. Valiviita. & Wandj 120061) . It is the same 
shape produced by the stand ard model of inflation 
dFalk. R angaraian. & Srednicki 1993) but in this case the am- 
plitude can be arbitrary. The potential bispectrum takes the simple 
form 

S,(P,,P 2 ,P 3 ) = 2A 2 /nl[(p 1 P2)"^ 4 + (PtPi)'"-- 4 + (P2P 3 )" S - 4 ] ,(19) 

and it is maximized for squeezed configurations. The quantity A is 
the spectral amplitude of the potential (given by <x 8 ), while rc s is the 
spectral slope. 

3.2 Equilateral shape 

This shape is a consequence of the inflaton Lagrangian being 
non-standard, and containing higher-order derivatives of the field 
(Alishahiha, Silv erstein. & Ton dl2004t lArkani-Hamed et al.ll2004 
iLi. Wa ng. & Wand l2008l) . The resulting bispectrum is maximized 
for equilateral configu rations. A template for the equilateral bispec- 
trum can be found in lCreminelli e t al. (2007), however the expres- 
sion is rather cumbersome and I did not report it here. The same 
applies to the following shapes. 

3.3 Enfolded shape 

The enfolded shape results from primordial non-Gaussianity be- 
ing evalu ated without the r egular Bunch-Dav i es va cuum hy- 
pothesis dChen et all 120071: lliolman & Tollevl l2008b . In this 
case the bispectrum is maximized for squashed configura- 
tions. A template for such a bispectrum is reported in 
iMeerburg. van der Schaar. & Cor asaniti (20091). 
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Figure 1. The shape of the function uim, quantifying the non-Gaussian con- 
tribution to the angular momentum variance growth. Different line styles 
and colors refer to different non-Gaussian bispectrum shapes, as labeled. 
The non-Gaussianity induced by the matter bounce is independent on /nl, 
while in all other cases /nl = 1 has been assumed. 

3.4 Orthogonal shape 

This shape is defined as being orthogonal (with respect to a suit- 
ably defined scalar product) to both the local and equilateral 
forms. The resulting bispectrum is maximized for both equilat- 
eral and squashed configurations, and a template can be found in 
ISenatore. Smith. & Zaldarriagal (120100 . 

3.5 Matter bounce shape 

This configuration is the consequence of a model universe without 
inflation, but with a scale factor that b ounces in a non-singular way 
dBrandenbergerl l2009; Cai et al 1 120091) . The matter bounce leads to 
a scale-invariant spectrum of density fluctuations, and to a bis- 
pectrum whose shape is similar to the local shape. Being non- 
inflationary in origin, the non-Gaussianity induced by a matter 
bounce model has no dependence on / NL . It can instead be shown 
that the matter bounce bispectrum is comparable to the local one 
with a fixed /nl = -35/8.1 considered explicitly the matter bounce 
because it is in principle possible that a weighted integral of the bis- 
pectrum, such as the one in Eq. l !17t , would magnify its differences 
with respect to the local model. As I show below this is actually not 
the case. 



4 RESULTS 

In Figure [T]l show oj m as a function of the mass scale for the five 
non-Gaussian cosmologies considered in this letter. In all cases, 
except for the matter bounce, I selected / NL = 1 in order to purely 
highlight the effect of the bispectrum shape. It is however easy to 
see that cj m is simply proportional to / NL . All curves decrease with 
increasing mass, implying that larger matter patches acquire lower 
amounts of angular momentum than smaller ones. This behavior 
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Figure 2. The non-Gaussian contribution to the spin variance growth nor- 
malized by the linear contribution, evaluated at the collapse time for an 
overdensity of a given mass. Line types and colors are the same as in Figure 
[T] and /nl = 1 has been assumed for all models except the matter bounce 
(which is /NL-independent). 

is similar to the mass dependence of the linear term given in Eq. 
{B). Also, curves referring to different models are rather similar 
in shape. Besides the matter bounce model, the local model is the 
one having the largest effect, while the orthogonal model has the 
lowest. This is different from the behavior of, e.g., the halo bias, 
for which the equilateral model displays the smallest effect. The 
function a> M for the matter bounce is virtually identical to that for 
the local model when assuming / NL = -35/8, in agreement with 
the previous discussion. 

As can be seen from the structure of Eqs. (1151 and 1 1161 . the 
relative importance of the non-Gaussian contribution with respect 
to the linear one is 



v , , _ „ </l(T)-/ 2 (T)> -dE(T)/dT U> M 



<IIJi(t)IP> 



dD(T)/dT d 



(20) 



The acquisition of angular momentum by protostructures happens 
at high redshift where, under the assumption of flat spatial geom- 
etry, the Universe is well approximated by an Einstein-de Sitter 
model. Corrections due to the presence of a cosmological constant 
can be considered to be negligible. If this is the case, then 



r M (T) = -^T- 2 ^ = -H (3 r) 2 ^ 



(21) 



Note that the negative sign cancels with the negative sign in the 
definition of a> M , so that the non-Gaussian contribution is positive 
(increases the spin growth) for a model with positive skewness and 
negative otherwise. This has already been noticed by CT. The pre- 
vious equation also shows that the non-Gaussian contribution to 
spin acquisition grows faster than the linear one. 

If the matter patch under consideration is an overdense re- 
gion, it is reasonable to assume that the spin growth induced by 
tidal torques occurs until the overdensity detaches from the over- 
all expansion of the Universe and collapses into a bound structure. 
This moment t, can be naively identified as D{T,)a M = 1- For an 



Einstein-de Sitter cosmological model this implies tJ = cr M , and 
thus 



T M (r,) = -H^ 
i en. 



(22) 



In Figure|2]l show the mass dependence of the function T m (t.) 
for the various non-Gaussian cosmologies that have been consid- 
ered in this work. For models with inflationary non-Gaussianity 
I assumed /nl = 1. As can be seen the mass dependence is in 
all cases relatively weak. In the local and matter bounce mod- 
els T m (t») is basically unchanged for masses ranging between 
M = 10 8 /r'M Q and M = \Q ls hr l M Q . For the equilateral and en- 
folded models Y m (t») increases by ~ 50% over the same interval, 
while for the orthogonal case it decreases by a factor of ~ 3. Hence, 
despite the fact that the linear and non-Gaussian contributions to the 
spin growth both decrease with mass, their relative importance re- 
mains relatively unchanged. The only exception is represented by 
the orthogonal model. 

Next, I selected a reference mass scale of M = W lo h' l M B 
and computed the dependence of T m (t») on /nl, shown in Figure 
[3] As previously mentioned, this dependence is always linear, how- 
ever the Figure is important in order to understand for what value 
of /nl a certain non-Gaussian model provides a given contribution 
to the total angular momentum variance. The matter bounce non- 
Gaussianity is independent of /nl, hence its contribution is always 
at the percent (negative) level compared to the linear one. As for the 
other models, in order for the non-Gaussian contribution to be com- 
parable to the linear one, primordial non-Gaussianity would need 
to be at the unrealistic level of /nl ~ 400 for the local case, and 
substantially larger than that for other bispectrum shapes. 

Figure [3] allows one to determine the non-Gaussian contri- 
bution, in units of the linear contrib ution, for the current bounds 
on /vjl- Constraints from the CMB dKomatsu et alj[20llh imply 
- 1 3 < / NL < 96 at 95 % Confidence Level (CL) for the local shapejj, 
meaning that the non-Gaussian contribution can be at most ~ 19% 
of the linear one. The same CMB data constrain -278 < /nl < 346 
for the equilateral shape, implying a ~ 18% relative importance. 
The tighter constraints on the le vel of non-Gaus sianity for the en- 
folded shape come from the LSS (Xia et al.l201 Th . corresponding to 
-16 < /nl < 465 at 2<x confidence level. This means that the non- 
Gaussian contribution is at most ~ 16% of the linear one. Finally , 
for the orthogonal shape the CMB data by iKomatsu et alj feOl ll) 
bear -533 < /nl < 8, implying a maximum ~ 10% (negative) rela- 
tive strength. These numbers can be appreciated also by looking at 
the positions of the filled circles in Figure [3] For comparison, the 
black dotted line shows the upper limit to the non-Gaussian con- 
tribution found by CT, after assuming a log-normal distribution for 
the primordial gravitational potential and after rescaling it to the 
scale M = 10 l0 h~ , M o . I stress the fact that, while CT adopted a 
Gaussian window function and h = 0.5, 1 used a real-space top-hat 
filter and h = 0.704. Moreover, while CT calibrated the level of pri- 
mordial non-Gaussianity using a value Sr = 4 for the skewness of 
the matter density field on a scale R = 8/r'Mpc, I adopted, conser- 
vatively, S R = 0. 1. This value results from the large-scale skewness 
per unit /nl for local non-Gaussianity (~ 10~ 3 , e.g., Figure 1 of 
iFedeli et al]|201 ll) multiplied by the most recent upper limit on the 
level of non-Gaussianity for the same shape (/ N l ~ 100). 



1 I converted all constraints mentioned her e to the LSS c onven tion, which 
has been adopted throughout. See IFedeli & Moscardini (2010) and refer- 
ences therein for a discussion. 



Angular momentum growth 5 




0.5 1 1.5 2 2.5 3 3.5 4 

L °g( f NL) 

Figure 3. The non-Gaussian contribution to the spin variance growth nor- 
malized by the linear contribution, evaluated at the collapse time for an 
overdensity of mass M = 10 10 /! _1 M o . Line types and colors are the same 
as in Figure [T] and results are shown as a function of /nl- The black solid 
line highlights the locus where the non-Gaussian contribution is identical 
to the linear one. Filled circles on each curve represent the maximum |/nlI 
values allowed by current constraints from CMB and LSS. The black dotted 
line shows the upper limit to the non-Gaussian spin contribution found by 
CT assuming a log-normal model for the primordial gravitational potential, 
after rescaling CT's result as detailed in the text. 

5 CONCLUSIONS 

I reconsidered the impact of primordial non-Gaussianity on the 
acquisition of angular momentum by CDM protostructures. Non- 
Gaussian initial conditions provide a next-to-linear correction to 
the spin growth that is absent when density fluctuations are nor- 
mally distributed. Previous results, obtained by CT after assuming 
a log-normal primordial gravitational potential, resulted in a contri- 
bution to the spin variance of ~ 10% with respect to the linear one. 
This value holds for a scale M = 10 10 /i _1 M G (with the cosmology 
of this work) and is based on a matter skewness of S r ~ 0.1, as 
deduced in the previous Section. Other models turned out to give 
a very large quasi-linear effect, suggesting that Lagrangian pertur- 
bation theory could not be successfully applied in those cases. I 
found that for several current models of non-Gaussian initial con- 
ditions, the contribution to the galactic spin variance during the 
mildly non-linear regime is similar to what predicted assuming a 
log-normal primordial gravitational potential. Considering the up- 
per limits to the current constraints on the level of inflationary pri- 
mordial non-Gaussianity returns a next-to-linear contribution at the 
level of ~ 10 - 20%. These results imply that the spin growth 
induced by inflationary non-Gaussianity seems to be generically 
tractable via perturbation theory. 

CT also demonstrated that higher-order contributions in the 
case of Gaussian density fluctuations provide a correction to the 
angular momentum variance equal to ~ 60% of the linear term. 
This means that the next-to-linear non-Gaussian contribution es- 
timated here has a significant impact on the spin acquisition by 
protostructures. Such an impact could potentially be even larger, 



because higher-order non-Gaussian corrections, that have not been 
considered here, depend on the trispectrum of the Zel'dovich poten- 
tial, and hence also react to primordial non-Gaussianity. The results 
presented in this letter motivate the study of these higher-order con- 
tributions, and show how it is possible to consistently describe the 
dynamics of protostructures based on largely general cosmological 
initial conditions. 

ACKNOWLEDGEMENTS 

I thank the University of Florida for support through the Theoret- 
ical Astrophysics Fellowship. I credit L. Moscardini for insight- 
ful comments on the manuscript and I am deeply indebted to the 
anonymous referee for help in substantially improving this work. 



REFERENCES 

Alishahiha, M., Silverstein, E., & Tong, D. 2004, Phys. Rev. D, 70, 123505 
Arkani-Hamed, N., Creminelli, P., Mukohyama, S., & Zaldarriaga, M. 

2004, Journal of Cosmology and Astro-Particle Physics, 4, 1 
Babich, D., Creminelli, P., & Zaldarriaga, M. 2004, Journal of Cosmology 

and Astro-Particle Physics, 8, 9 
Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 

15 

Bartolo, N, Komatsu, E., Matarrese, S., & Riotto, A. 2004, Phys. Rep., 
402, 103 

Bernardeau, F. & Uzan, J. 2002, Phys. Rev. D, 66, 103506 
Brandenberger, R. 2009, Phys. Rev. D, 80, 043516 

Cai, Y.-F, Xue, W., Brandenberger, R., & Zhang, X. 2009, J. Cosmology 

Astropart. Phys., 5, 11 
Catelan, P. 1995, MNRAS, 276, 115 
Catelan, P. & Theuns, T. 1996a, MNRAS, 282, 436 
Catelan, P. & Theuns, T. 1996b, MNRAS, 282, 455 
Catelan, P. & Theuns, T. 1997, MNRAS, 292, 225 
Chen, X. 2010, Advances in Astronomy, 2010 

Chen, X., Huang, M., Kachru, S., & Shiu, G. 2007, Journal of Cosmology 

and Astro-Particle Physics, 1, 2 
Creminelli, P., Senatore, L., Zaldarriaga, M., & Tegmark, M. 2007, Journal 

of Cosmology and Astro-Particle Physics, 3, 5 
Doroshkevich, A. G. 1970, Astrophysics, 6, 320 
Falk, T., Rangarajan, R., & Srednicki, M. 1993, ApJ, 403, LI 
Fedeli, C, Carbone, C, Moscardini, L., & Cimatti, A. 2011, MNRAS, 

414, 1545 

Fedeli, C. & Moscardini, L. 2010, MNRAS, 405, 681 
Heavens, A. & Peacock, J. 1988, MNRAS, 232, 339 
Holman, R. & Tolley, A. J. 2008, Journal of Cosmology and Astro-Particle 
Physics, 5, 1 

Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 

Li, M., Wang, T., & Wang, Y. 2008, Journal of Cosmology and Astro- 
Particle Physics, 3, 28 

Meerburg, P. D., van der Schaar, J. P., & Corasaniti, S. P. 2009, Journal of 
Cosmology and Astro-Particle Physics, 5, 18 

Moscardini, L., Matarrese, S., Lucchin, R, & Messina, A. 1991, MNRAS, 
248, 424 

Peebles, P. J. E. 1969, ApJ, 155, 393 
Peebles, P. J. E. 1971, A&A, 11, 377 

Sasaki, M., Valiviita, J., & Wands, D. 2006, Phys. Rev. D, 74, 103003 
Schafer, B. M. 2009, International Journal of Modern Physics D, 18, 173 
Senatore, L., Smith, K. M., & Zaldarriaga, M. 2010, Journal of Cosmology 

and Astro-Particle Physics, 1, 28 
Shandarin, S. F. 1980, Astrophysics, 16, 439 
Sugiyama, N. 1995, ApJS, 100, 281 
White, S. D. M. 1984, ApJ, 286, 38 

Xia, J.-Q., Baccigalupi, C, Matarrese, S., Verde, L., & Viel, M. 2011, J. 

Cosmology Astropart. Phys., 8, 33 
Zel'dovich, Y. B. 1970, A&A, 5, 84 



